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Abstract. The synchrony of electric power systems is important in order to maintain stable 
electricity supply. Recently, the measure basin stability was introduced to quantify a node’s 
ability to recover its synchronization when perturbed. In this work, we focus on how basin 
stability depends on the coupling strength between nodes. We use the Chilean power grid 
as a case study. In general, basin stability goes from zero to one as coupling strength 
increases. However, this transition does not happen at the same value for different nodes. 
By understanding the transition for individual nodes, we can further characterize their role in 
the power-transmission dynamics. We find that nodes with an exceptionally large transition 
window also have a low community consistency. In other words, they are hard to classify 
to one community when applying a community detection algorithm. This also gives an 
efficient way to identify nodes with a long transition window (which is computationally time 
consuming). Einally, to corroborate these results, we present a stylized example network with 
prescribed community structures that captures the mentioned characteristics of basin stability 
transition and recreates our observations. 


PACS numbers: 64.60.aq, 84.40.Az, 89.70.-a, 89.75.Pb 
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1. Introduction 

Electric power systems are often modeled as oscillators on networks [[II El El . The rotational 
motion of power generators for alternating currents needs to be synchronized with a rated 
frequency of an electric power grid [|4l[5l[6l|71|8l|9l. Stable synchrony of electric power 
grids’ nodes is important to prevent cascading failures. When a node gets an external 
perturbation leading the node away from a synchronous state, then the node either returns 
to the synchronous state or escapes from its basin of attraction, typically to a different limit 
cycle. Basin stability is calculated as the fraction of the possible phase values a node can be 
perturbed to and still recover synchronicity [[^. It has become the standard way to quantify 
a power-grid node’s robustness to large point perturbations. The basin stability is determined 
by the magnitude of the perturbation and network characteristics of the node. It also depends 
on the coupling strength between the nodes in the dynamic system representing the electric 
power flow. For instance. Refs. |j71 [8l investigate topological characteristics of particularly 
stable or unstable nodes. The authors conclude that, in general, dead ends weaken the basin 
stability, while detours strengthen it. Such simple indirect estimates are helpful because basin 
stability is very computer intensive to calculate. 

These previous studies, however, study basin stability for one specific coupling 
strength [jbl El El. However, basin stability is sensitive to the coupling strength. In a 
power grid, the transmission capacity, which determines to the coupling strength, is decided 
based on the physical conditions such as the air temperature, the length of transmission 
line, and the amount of transmitting electricity. When the amount of electricity transmitted 
exceeds the transmission capacity, it may cause the breakdown of a large part of the system. 
Therefore, it is important to understand the functional dependence of basin stability on the 
coupling strength, not only its behavior at a fixed and intermediate coupling strength. We 
use the real electric power grid in Chile as our prime example, but also discuss how the 
results generalize. To be more specific, we investigate the network-structural factors [TTOlfTTII 
behind a node’s basin stability as a function of the coupling strength. Based on the Chilean 
case study, we argue that there is a connection between the coupling strength dependence 
and community structure—the way the network can be decomposed into communities that 
are densely connected within and sparsely connected between each other [Oil . Other 
authors have noticed that those communities easily synchronize internally. In a similar spirit, 
we investigate how the membership strength of a node within a community decomposition 
(its community consistency) of the network relates to the coupling strength dependence, 
specifically the width in parameter space of the basin stability’s transition from zero to one. 
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2.1. The dynamical model of electricity transmission 

The synchronization dynamics between power-grid nodes is commonly modeled as the set of 
nonlinear oscillators using a Kuramoto-type model [[Ill2lllll5ll6l|7l[8l|9l|: 



( 1 ) 


where Pi is the net power generation (positive) or consumption (negative) at node z; cr, is a 
dissipation constant; K is the coupling strength between nodes i and j; the adjacency matrix 
Oij = 1 if there is an edge between node i and j, and a,y = 0 otherwise; 9i is the phase of 
node z; a>i = 6i is the angular velocity of node i. A power grid is considered to be stable and 
synchronized when cu, vanishes for all of the nodes so that the system maintains the desirable 
constant frequency. For numerical integration, we use the Runge-Kutta method [|T^ with the 
convergence criterion (Til of the time derivative of angular frequency zb < 5 x 10~^ and the 
angular frequency < 5 x 10“^ based on the actual fluctuations, i.e., we consider the initial 
phase and angular velocity (0„ uj,) of z as belonging to the stable basin if the system eventually 
converges to tu < 5x 10“^ and to < 5x 10“^ after the z is perturbed. We perturb z by initializing 
tOi and di to random values in the intervals [-100,100] and {-n, ;7r]. The other nodes have their 
values of the synchronous state. Following Refs. [I^|7l[8l, we set or, =0.1 for all z. 

2.2. The Chilean power grid 

We construct the Chilean power grid, our prime example, from the operational records 
published by the National Energy Committee of Chile liTTlI and the data regarding the power- 
grid transmission lines published by Centro de Despacho Economico de Carga del Sistema 
Interconectado Central (CDEC-SIC)—the major electric power company in Chile [fTSl . The 
Chilean power grid consists of power plants and substations as nodes, and the transmission 
lines between them as edges IfTOl fTTII . A power plant generates electricity and a substation 
distributes it to the final consumers. We aggregate the power grid using Zhukov’s aggregation 
scheme IfIPl . resulting in nodes being either net generators (F, > 0 in Eq. Q) or net consumers 
{Pi < 0) based on the real data [fl^ . After determining the sign of F, values, we first set 
Pi = -1 for Pi < 0 and F, = Fpos > 0 so that YiiPi = 0 (Kirchhoff’s circuit laws). The 
inactive nodes such as towers or net zero traders (F, = 0 in Eq. ([T])) do not intervene with the 
perturbation propagation and are removed. We apply the so-called Kron reduction [|^ 1211 
procedure. It is a node elimination technique widely used in electrical circuit design. Eor a 
given network configuration, Kron reduction removes a node and adds new edges between 
neighbors of the eliminated node at each step [l22l l23l . Consequently, as the more nodes 
are Kron-reduced, the denser (in terms of mean degree) the network becomes. However, 
even though the network structure is changed, the new edges are placed so that the network 
maintains the relationship between nodes in terms of electricity flow. Eor the Chilean power 
grid, we eliminate 46 inactive nodes using Kron reduction. As a result, the final network 
consists of 291 net consumers, 129 net generators, and 573 edges. 
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Figure 1. Color-map and distributions of basin stability values of nodes in the Chilean power 
grid for different coupling strength K - 1.2710 (a), 1.2715 (b), 1.2720 (c), and 1.2725 (d). 


2.3. Basin stability transition 


Basin stability is a stability measure for a network [IH Ul [8]|. For eaeh node, an initial 
phase value is seleeted uniformly at random. Then we run the dynamies until eonvergenee 
(aeeording to the eriterion presented in See. ^A_). The fraetion of initial phase values that 
leads to eonvergenee defines the basin stability. We sample 100 points uniformly at random 
from the intervals -100 < m < 100 and -n < 6 < n. Assuming a uniformly random sampling 
is done for the sake of simplieity (following Refs. [I6l|71|8l). Note that the basin stability of a 
node eould be eompared to other nodes in the same network, but not neeessarily to nodes in 
other networks. 

For the so-ealled infinite bus-bar model (basieally a mean-field approximation of the 
problem we study), the basin stability inereases gradually with the transmission eapaeity and 
suddenly beeomes unity (stable) [[61. This jump appears when the transmission eapaeity is 
large enough for all perturbations or initial shoeks from the entire phase spaee to be absorbed 
into the system. However, in the more realistie multi-node model, the basin stability and its 
transition form are different for different nodes. So for a eomplete pieture, every node needs 
to be analyzed separately. 
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Figure 2. A schematic diagram of basin stability transition as a function of coupling strength 
on a power grid. 


We show the results for the distribution of basin stabilities of the Chilean power grid’s 
nodes as funetions of the eoupling strength ^ in fig. For low K, all nodes have zero basin 
stability. As K inereases, a few nodes start to reaeh a non-zero basin stability (fig. 

When the eoupling strength reaehes a eertain threshold, almost all the nodes suddenly turn 
into a stable state. Note that the basin stability shows a nonlinear transition as a funetion of 
eoupling strength. The dramatieally different distribution of basin stabilities for the different 
K reinforees the importanee of eonsidering not only a eertain basin stability at a ^ value but 
also the transition of basin stability eorresponding to different K values. In See.|^ we address 
the transition of basin stability and speeifieally introduee basin stability transition window. 

3. Results 

3.1. Basin stability transition window 

Basin stability ehanges, as mentioned, with the eoupling strength. Figurej^shows a sehematie 
diagram of the basin stability transition pattern that we eall transition window, for two nodes 
of a hypothetieal network. At ^ = Ki^, node 1 has larger basin stability than node 2. On the 
other hand, when the eoupling strength inereases io K = Ki, node 2 beeomes more stable than 
node 1. These transitions of basin stability of nodes 1 and 2 are represented as a funetion of 
K on the transition window (the shaded reetangles of fig.[^. 

To loeate the transition windows of the nodes, we use the biseetion method ^2M . We 
start by finding the lower bound of the transition window ^low, whieh is the minimum value 
of K that makes the basin stability larger than zero. We do this by traeking the basin stability 
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Figure 3. Basin stability transition of node A, B, and C (see fig.j^for the actual locations of 
the nodes) (a) and average transition windows (b). In (a), the solid gray lines are the guide to 
the eyes. The dots represent the average basin stability of every ten points with the resolution 
0.1 from = 0 to 19.9 (each dot located at K corresponds to the averaged values of basin 
stability for K - 0.5, K - 0.4, ..K + 0.4), where error bars corresponding to the standard 
deviation for each of the uniform bins with the size 1.0. Nodes A, B, and C reach the complete 
basin stability (i.e. unity) at = 1.3, 13.1, and 14.7, respectively. In fig. [^b), on the left panel, 
the fraction of nodes belonging to each logarithmic bin is shown for each bin. On the right 
panel, the average AK and are shown for each bin, where the error bars show the standard 
deviation values of /flow (on the left) and /Thigh (on the right). The inset shows the AK and /Tmid 
values in a transition window (the filled rectangle) for node B. 


from the ^min = ^min,init hnd K^ax = ^max.init, and measure the basin stability of a node at the 
midpoint i^midpoint = (^min + ^max)/2. When the basin stability at i^midpoint is non-zero, i^max is 
set to the eurrent ^midpoint value and the new ^midpoint value is reealeulated from the new ^max 
value. Similarly, when the basin stability at ^midpoint is zero, the new ^min is set to the eurrent 
^midpoint- This iteration eontinues until < ^threshold and K\o^ is set to the final values 

of (.fiTmin + ^max)/2. The Upper bound .fiThigh of the transition window (the maximum value of 
K that makes the basin stability smaller than unity) is also determined in the same fashion. 
As a result, we obtain the transition width lAK = ^high - ^low that represents the width of the 
transition window. We use i^min.init = 0, i^max,init = 20, and i^threshoid = 0.01. 

ls.K can vary much between the nodes. For instance, node A (one of the nodes with the 
narrowest lS.K values) of the Chilean power grid shown in fig.|^a) undergoes the sudden basin 
stability jump at ^ - 1.3 and reaches the completely stable state at a smaller K value than 
other nodes such as nodes B and C as comparison. The distribution of lAK in the Chilean power 
grid is highly right-skewed (with a mean of 0.154 and median of 0.022). Figure [^b) shows 
various statistics of the transition window, where we divide the lAK values into five different 
logarithmic bins and present the fraction of nodes belonging to each bin (the left panel of 
fig. [^b)) and the range of average [^low, ^high] for each bin (the right panel of fig. [^b)). The 
vast majority (about 96%) of nodes have very small values of lAK <0.1 (about 35% of nodes 
even have ls.K < 0.01). However, there are a few nodes with relatively very large lAK values, 
e.g., the average range [^low, ^high] for the bin corresponding to the maximum K range (10% 
10^) is more than 13 times as large as that for the bin corresponding to the minimum K range 
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(10~^-10~^). In general, the wider of the transition window a node is, the larger K value in 
the middle ((^low + ^high)/2, heneeforth denoted as ^mid) the node has. 

It is interesting to note that the transition curves themselves seem to have different shapes 
for different nodes, even for the nodes with similar tsK values. For example, l^K of nodes B 
and C have the similar tsK values but the former (latter) node shows roughly concave (convex) 
curves at the transition. The different range of tsK and the shape can be interpreted in several 
ways. If a node has small IsK, the location of is also small according to our observations, 
so the synchrony of the node tends to be stable even for small values of coupling strength. In 
other words, the node stays functional and is hardly affected by the accidental performance 
drop of the transmission line connected to it. However, at the same time, once the basin 
stability value started to decrease, the node suddenly loses the stability. Consequently, it 
makes the detection of system failure difficult. When the transition of the basin stability 
is gradual, on the other hand, a system disturbance can be detected by a drop in the basin 
stability. In summary, narrow transition windows could be good for keeping stability, while 
the wider window are better for an early detection of system failure. 

3.2. Chilean power-grid analysis 

One drawback of basin stability as a robustness measure of nodes in power-grids, is that 
it is computationally demanding. Therefore, it is desirable to find less computationally 
restrictive, indirect measures to estimate the basin stability. Some previous studies proposed 
such topological indicators for basin stability—Menck et al. [|3 predict that dead-end nodes 
have small basin stability, while Schultz et al. [[81 predict that detour nodes have large basin 
stability. 

The basin stability estimators in Refs. 0 and 0 assume a specific .fif-value. If K is 
not precisely known or varying, it can only identify extreme cases (very stable or unstable 
nodes). The width of transition window is thus a more appropriate measure to capture the 
role of a node in the transmission dynamics. To understand why some nodes have larger ls.K 
than others, we try to find the most appropriate explanatory measure. This is far from trivial. 
Indeed phase synchronization on a network is a complex consequence of the propagation of 
phase difference and recovery due to the interactions between nodes 0. 

It is known that the mesoscopic properties of networks, such as community structure, 
play a significant role in the synchronization on networks [[251 [26l [27l . A community is a 
subnetwork, that is more strongly connected within than to the rest of the network ffT^ [T3l . 
Different community detection methods could divide the same network in different ways, 
but there would typically be pairs of nodes that always are classified as belonging to the 
same community. Some community detection methods are non-deterministic and could also 
give a different community decomposition between different runs of the same algorithm. 
One could therefore define nodes often classified to the same community as having high 
community consistency with respect to this algorithm. We use GenLouvain [[28ll . a variant 
of conventional Louvain community detection algorithm Il29ll known for its computational 
scalability, to identify communities in networks. We can also control the resolution parameter 
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Figure 4. Four different community detection results are shown for the same resolution 
parameter y - 0.0323, where the regions of nodes with inconsistent community assignments 
are magnihed. Different colors represent different communities detected by the GenLouvain 
method. 

in the modularity function used as the objective function to maximize, to detect communities 
with different scales [|^ 1311 133 • We define community consistency of individual nodes as 
the nodes’ degree of certainty of the assigned community. Previously, Refs. [fT?l[T5]| calculate 
community consistency (or “consensus clustering”) among assigned communities in multiple 
runs of each stochastic community detection algorithm, where they focus on the community 
consistency of each algorithm. In this paper, however, we measure community consistency of 
individual nodes to relate it to the nodes’ basin stability transition width ls.K. 

Figure shows four different community detection results of the GenLouvain method 
on the Chilean power grid, where the we use the same resolution parameter OOl [311 l32l — 
y = 0.0323. We chose this value to get around three communities which matches our visual 
impression of the network. As hinted in fig.|^ most nodes are assigned to the same community 
for all iterations. We call such nodes community consistent. Some nodes, however, are 
assigned to different communities in different runs (the magnified region in fig. |^. In terms 
of power-grid dynamics, such nodes could, be believe, be influenced by perturbations from 
different directions (communities). 

To measure the community consistency of individual nodes, we first observe a node i 
is perfectly community consistent if, and only if, it gets always classified into a cluster with 
the same other nodes. Let 0,, be the fractions of runs of the community detection algorithm 
when i and j are classified into the same community, then i is perfectly community consistent 
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Community consistency 


Figure 5. AK and community consistency. This is basically a scatter plot but darker points 
represent more points, drawn with transparency. 


Table 1. Pearson correlation coefficient r of AK versus community consistency (<!)), degree 
(k), clustering coefficient (C), and current flow betweenness (F) centrality. 



O 

k 

C 

F 

r 

-0.581 

0.033 

-0.054 

0.072 

p-value 

< 10-3 

0.500 

0.266 

0.139 


if and only if is either 1 or 0 for all j. Sinee (pij e [0,1], we ean obtain a metries for 
community consistency by measuring the average distance to (pjj = 1 /2. But rather than the 
linear distance (|0,y - 1 /2|), we sum the square distance. First, this conforms our measure to 
standard measures of deviation or spread like the root-mean-square, variance, and radius of 
gyration. Second, the non-linearity of the parabola accentuates values close to 0 or 1 and tones 
down the mid-interval values. This turns out to be practical since most 0-values are close to 0 
or 1. Finally, we multiply a factor 4 to the outlined measure to get a value in the unit interval. 
In summary, the community consistency of i is given by 

< 2 ) 

j*i 

The correlation between community consistency and ls.K of basin stability transition 
window of nodes is shown in fig. and Table where the larger values of community 
consistency a node has, the narrower lAK the node has. In other words, if a node’s community 
consistency is weak, so that the node is assigned to different communities for each different 
realization, the node tends to have a wide basin stability transition window. Figurej^illustrates 
the nodes where lAK and community consistency is colored, which provides a visual evidence 
of our result of the correlation between lAK and community consistency. It seems that the 
basin stability transition occurs at some characteristic values of coupling strength for different 
communities of nodes as a unit, and the transition threshold is somehow smeared out for those 
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Figure 6. hKj (the left panel) and community consistency (the right panel) of the 
Chilean power grid, where A/Tmax is the maximum IsK value among the nodes. The insets 
show the area of interest where the nodes with large AK and small community consistency. 


nodes with weak or inconsistent community membership, which makes IsK large. There are 
several nodes in the northern part of the power grid with lower community consistency, but 
we believe that the effect of the nodes is not strong enough to affect basin stability transition 
window. Moreover, compared to the actual numerical integration for calculating the basin 
stability, the calculating community consistency is much faster (at least for reasonably fast 
community detection algorithms). 

To check if there is any other (possibly simpler) network measures that can predict 
lS.K, we measure other network metrics. The Pearson correlation coefficients for lS.K versus 
degree ifTOl [TT]| and versus community consistency are about 0.033 (p-value around 0.5) 
and -0.581 (p-value less than 10“^), respectively. For the correlation coefficient values for 
other representative network centralities and community consistency are shown in Table 
indicating that only the mesoscopic measure of community consistency is significantly 
correlated with lS.K, in contrast to other conventional measures such as degree, clustering 
coefficient ttni, and current flow betweenness centralities (supposedly a more relevant type of 
betweenness in our power-grid case than the ordinary one) [fTTl . The two former centralities 
are microscopic or local, while the latter is macroscopic as it deals with all of the possible 
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Figure 7. Correlation between /S.K and degree (a), clustering coefficient (b), and current flow 
betweenness (c) centralities. These are basically scatter plots but darker points represent more 
points, drawn with transparency. 


pathways for an entire network. We tried other measures than the ones presented in Table 
but we eould not find anything statistieally signifieantly eorrelated with ls.K. 

3.3. Example network analysis 

In order to verify our results from the Chilean power grid that the basin stability transition 
is elosely related to the eommunity membership (figs, [^a) and (b)) we eonstruet a simple 
example network with preseribed eommunity struetures depieted in fig. 8(o). The example 
network eonsists of 18 nodes and most nodes are separated into two eommunities eaeh of 
whieh has six fully eonneeted nodes and a single attaehed dead-end node, exeept for four 
nodes branehed from the bridge between the eommunities. The square nodes are assigned 
as eonsumers (P, = -1 in Eq. Q) and the eireular nodes are produeers (E, = 1 in Eq. Q). 
The example network symbolizes a network strueture with two well-defined eommunities and 
outlier nodes on an interfaee braneh. The interfaee braneh represents a eonneeted subgraph 
without elear eommunity membership, just like the nodes E and G in the Chilean power grid 
(shown in fig.j^a)). We measure basin stability transition window of the nodes in the example 
network and show representative eases in fig. [^d). Most eharaeteristies of the transition 
patterns observed from various nodes in the Chilean power-grid nodes (fig.j^b)) are eaptured 
by some representative nodes in the example network (fig.j^d)), e.g., eoneave versus eonvex 
transitions and the stairease pattern. 

To be more speeifie, the basin stability transition eurve beeomes eonvex or shows the 
stairease pattern, when the node is loeated far from a preseribed eommunity and near to the 
braneh from the bridge (thus with a weak eommunity membership), in both Chilean power 
grid and our example network. Sueh outlier nodes seem to have a tendeney of reaehing a stable 
state (the basin stability equals to unity) at larger values of K eompared to the nodes with 
stronger eommunity membership. In other words, our example network ean be eonsidered 
as the simplified version of the Chilean power grid at least in terms of the basin stability 
transition. In partieular, we find some speeial transition patterns in both networks. Node G in 
the Chilean power grid (figs.j^a) and (b)) shows a stability at a smaller value of K eompared 
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Figure 8. Basin stability transition windows of some nodes in the Chilean power grid, (a) and 

(b) , and the example network, (c) and (d). Same types of nodes in terms of basin stability 
transition patterns for the Chilean power grid and our example network, according to our 
judgment, are denoted with the same color, both for the nodes in (a) and (b), and the curves in 

(c) and (d). Numbers on the nodes are used for the identification in the corresponding basin 
stability transition plot, for each network separately. 


to Other nodes in the bridge. When a node is loeated near to the terminal node, the transition 
eurve ehanges its shape form from eoneave down to eoneave up. However, the aetual terminal 
node G in the Chilean power grid does not follow this trend and reaehes medium level of basin 
stability at relatively small K and maintains the plateau, until the basin stability suddenly 
reaches unity. Node 6 in our example network (fig. [^d)) shows a similar pattern. Those 
nodes are commonly located at the terminus of the interface branch, and we speculate that 
this pattern is related to that the nodes cannot spread external perturbation from the neighbor 
nodes to the rest of the community. There are also some differences between the Chilean 
network and the example network, for example the spikes from the plateau to the maximum 
basin stability (fig.j^d)). For example, nodes 1 and 6 in the example network (figs.j^c) and 
(d)) has maximum basin stability at ^ « 5 and 10 (node 1) and ^ « 13 and 17 (node 6), and 
they become unstable again. We find this phenomenon for the small example network too (not 
shown). Why this happens is beyond the scope of this paper, but not surprising in the light of 
the instabilities in synchronization on networks [13| . 
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We have studied the basin stability transition throughout the eoupling strength parameter 
spaee, foeusing on the basin stability transition window as a new metries for eharaeterizing 
the eontribution of a node to the stability of a power grid. As previous works focus on the 
stability for specific coupling strength values, our approach complements these. While a 
narrow transition window implies a sudden change in the stability, a wide transition window 
can provide an early signal of danger when the coupling strength is gradually weakened. By 
comparing the mesoscopic network property of community consistency with the transition 
window, we found that the former is a good predictor of the latter (or vice versa), signifying 
the importance of community structures on the synchronization dynamics again. 

On a practical side, such community-consistency based predictions provide a proxy 
of the actual time-consuming simulations of calculating basin stability, especially for very 
large systems. Furthermore, once we assign the communities and the strength of nodes’ 
memberships by community consistency, we are able to analyze the dynamics in the unit 
of such communities instead of the entire node set. We would also like to emphasize that 
the network example with given community structures and the bridge effectively captures the 
basin stability transition properties observed in the Chilean power grid. 

Our approach can be improved by considering more realistic situations such as assigning 
different parameters of real fluctuations of power input and output, different edge weights 
based on admittance of the transmission line, etc. We see our work as a starting point, and 
anticipate more studies in the future. Finally, we would like to emphasize that our community- 
consistency measure is not limited to the power-grids, but could be used in all kinds of 
community detection problems. 
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